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Abstract: The production of doubly charged Higgs bosons (-ff ±± ) at the CERN LHC 
can give rise to distinctive multi-lepton signatures. The discovery potential of ff ±=t can 
be optimized by considering a search strategy which is sensitive to both of the dominant 
production mechanisms, qq — > H ++ H and qq' — > H^H^ . We compare the discovery 
potential for the signatures of exactly four leptons and at least three leptons in the final 
state, using the same set of cuts. We have carried out fast detector simulations at the LHC 
for both signal and backgrounds for a wide range of values of the charged Higgs mass. We 
find that the use of the latter channel can substantially improve the detection prospects of 
the doubly charged Higgs boson at the LHC. 
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1. Introduction 



The firm evidence that neutrinos oscillate and possess small masses below the eV scale Q 
necessitates physics beyond the Standard Model (SM), which could manifest itself at the 
CERN Large Hadron Collider (LHC) and/or in low energy experiments which search for 
lepton flavour violation (LFV) Q. Consequently, models of neutrino mass generation 
which can be probed at present and forthcoming experiments are of great phenomenological 
interest. 

Neutrinos may obtain masses via the vacuum expectation value (vev) of a neutral Higgs 
boson in an isospin triplet representation [Q, |], ||, [6|, [7J. A particularly simple implementa- 
tion of this mechanism of neutrino mass generation is the "Higgs Triplet Model" (HTM) 
in which the SM Lagrangian is augmented solely by an SU{2) triplet of scalar particles 
with hypercharge Y = 2 ||, ^, In the HTM, neutrinos acquire Majorana masses given 
by the product of a triplet Yukawa coupling (hij) and a triplet vev (v a). Consequently, 
there is a direct connection between and the neutrino mass matrix, which gives rise 
to phenomenological predictions for processes which depend on h tj @, g, [fi], 0, ||, ||, 



14, 15, 16, 17, 18, [L9|. A distinctive signal of the HTM would be the observation of a 



doubly charged Higgs boson (ff ±± ), whose mass (M H ±±) may be of the order of the elec- 
troweak scale. Such particles can be produced with sizeable rates at hadron colliders in 
the processes qq -)• H++H— pi, P, [2l M, p|] and qq 7 H^H^ 



2C, p^, 26 1, where 



is a singly charged Higgs boson in the same triplet representation. Direct searches for 
H ±ziz have been carried out at the Fermilab Tevatron, assuming the production channel 
qq — >• H ++ H and the leptonic decays ff ±± — > if if (i = e,^,r), and mass limits in the 
range M H ±± > 110- 150 GeV have been obtained p7|, M, M M . The CERN Large Hadron 
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Collider (LHC), using the above production mechanisms, will offer improved sensitivity to 



m h±± tea, papa 31, u 



The decay channels — > £~^£~^ and — > l^v are the dominant ones if v& ^ 10 -4 
GeV, and give rise to multi-lepton signatures. In the HTM, one expects ^ 10 -4 GeV 
if the triplet Yukawa coupling is larger than the smallest Yukawa coupling in the SM 
(i.e., the electron Yukawa coupling). One can define various multi-lepton signatures which 
originate from the production mechanisms qq — > H ++ H and qq' — > H^H^ . The four- 
lepton signature (4^) only receives a contribution from qq — > H ++ H , and the detection 
prospects in this channel at the LHC have been studied in ||, H ||, ||, [l| . Although 



this 4£ signature provides a very promising way to search for it is not necessarily the 

channel which offers the best sensitivity for a given integrated luminosity and mass M H ±± . 
Recently, attention has been given to the three- lepton channel |l4], which also has 
relatively small SM backgrounds. Importantly, the signature of three-leptons is sensitive 
to the production mechanism qq' — > H^H^ p6fl . The magnitudes of the cross sections 
of qq — > H ++ H and qq' — > H^H^ are comparable in a large parameter space of 
the HTM, because the scalar potential of the model gives M H ±± ~ M H ± (unless a specific 
scalar quartic coupling is taken to be fairly large) . In order to improve the sensitivity of the 
LHC to M H ±±, one can define two distinct signatures consisting of three leptons, in which 
two of the leptons have the same electric charge and are assumed to originate from ff ±=t 
(which is produced in both production mechanisms, qq — > H ++ H and qq' — > H^H^): 
i) the signature of "exactly three leptons" (3£), [15], and ii) the signature of "three or more 
leptons" (> 3£) p9| . Different sensitivity to M H ±± is expected in these two channels. 

Detection prospects at the LHC are best for t = e, [i, for which there are several 
exclusive three-lepton channels, e.g. eee, /i/i/i, e/i/i etc |l8| . However, in the region of high 
invariant mass for a pair of same-sign leptons (e.g., mu > 200 GeV) one expects similar 
detection efficiencies and SM backgrounds in these exclusive three-lepton channels with 
£ = e,fj,. Therefore, to optimise the sensitivity to M H ±± it is reasonable to define an 
inclusive signature in which and ^ are treated as identical particles, and this is the 
approach which is taken in [15]. The signature of 3£ (with £ = e,\x) is studied in [15] 
and it is shown that the sensitivity to M H ±± is significantly superior to that of the A£ 
channel. This very promising result can be further improved, because the > 3£ channel is 
expected to give even greater sensitivity to M H ±± than the 3£ channel. In this work we 
study the signature of > 3£ and compare its sensitivity to M H ±± with that obtained for the 
4£ channel. Our study is the first simulation of the signature > 3£ at LHC that includes 
both production mechanisms qq — > H ++ H and qq' — > H ±ziz H^. We show that the > 3£ 
channel is the optimum search strategy for H ±ziz . We note that the most recent search 
by the DO collaboration [29] uses the strategy of > 3£ in the context of a search for three 
muons (/i ± /i ± /i =F ), with the assumption that production of H ±±L is only by qq — > H ++ H . 

The paper is organized as follows. In section Q we will briefly review the HTM model. 
In section || we discuss the production of doubly charged Higgs bosons at hadronic colliders. 
In the same section we will describe our analysis setup and framework for simulations of 
backgrounds and signal processes. The results of our simulations are given in section 
Finally we will conclude with a summary of our results in section |5] 
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2. The Higgs Triplet Model 



a = :, . ^ . (2.2) 



The HTM model |, |, is an extension of the SM in which only the scalar sector is 
augmented with a Higgs triplet. It is a particularly simple model which contains a dou- 
bly charged scalar. The model has the following SU(2) <g> U(1)y gauge-invariant Yukawa 
interactions: 

£ 3 hi j ipf L Cia 2 Aipj L + h.c. , (2.1) 

where the triplet Yukawa couplings hij(i,j = e,[A,r) are complex and symmetric, C is 
the Dirac charge conjugation operator, 02 is a Pauli matrix, ViL = ( u iik)~L is a ^ e &- 
handed lepton doublet, and A is a 2 x 2 representation of the Y = 2 complex triplet fields 
(S++,S+,S°): 

( 5+/V2 5 ++ \ 
^ 5° -5+/V2J 

Note that the mass eigenstate i^ ±± is entirely composed of the triplet field (i/ ±=l= = <5 ±± ), 
while is predominantly 5 , with a small component of isospin doublet scalar (<!>). A 
non-zero Higgs triplet VEV, {5°} = v A /V2, gives rise to the following Majorana mass 
matrix for neutrinos: 

mij = 2h ij {5°) = V2h ijVA . (2.3) 

Realistic neutrino masses can be obtained with a perturbative hij provided that v a ^ 1 
eV. The presence of a non-zero va gives rise to p 7^ 1 at tree level, where p = / (M| cos 2 Ow)- 
Therefore ^ 1 GeV is necessary in order to comply with the measurement of p ~ 1. 
This simple expression of tree-level masses for the observed neutrinos is essentially the 
main motivation for studying the HTM. It provides a direct connection between and 
the neutrino mass matrix, which gives rise to phenomenological predictions for processes 
which depend on hij §, |, 0, 0, 0, 0, [b§ [jj], 0, [L7|, @, @- 



The mass matrix for three Dirac neutrinos is diagonalized by the PMNS (Pontecorvo- 
Maki-Nakagawa-Sakata) matrix Vpmns [0- F° r Majorana neutrinos (which is the case in 
HTM), two additional phases appear, and then the mixing matrix V becomes 

V = Vpmns x diag(l, e^/ 2 , e^ 2 ), (2.4) 



where 4>i an d 4>2 are referred to as the Majorana phases |g, 35] and — ir < 4>i,4>2 < One 
has the freedom to work in the basis in which the charged lepton mass matrix is diagonal, 
and then the neutrino mass matrix is diagonalized by Vpmns- Using Eq. ( ^^ ) one can 
write the couplings hij as follows || ||: 



_ rriij = 1 
V^A V2.VA 



FpMNsdiag(mi,m 2 e^m 3 e^ 2 )Vp T MN s . (2.5) 



Here mi, 7712 and 7713 are the absolute masses of the three neutrinos. Neutrino oscillation 
experiments are sensitive to mass-squared differences, Am|i(= m| — m\) and Am|i(= 
m| — m 2 ). Since the sign of Am^ is undetermined at present, distinct patterns for the 
neutrino mass hierarchy are possible. The case with Am^ > is referred to as normal 
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hierarchy (NH) where mi < 7712 < m.3, and the case with Am^ < is known as inverted 
hierarchy (IH) where m.3 < mi < m^- However, information on the mass rriQ of the lightest 
neutrino (either m\ or 7713) and the Major ana phases cannot be obtained from neutrino 
oscillation experiments. This is because the oscillation probabilities are independent of 
these parameters, not only in vacuum but also in matter. An attractive feature of the 
HTM is this simple relationship between the triplet Yukawa couplings and the parameters 
of the neutrino mass matrix (many of which are measurable) given in Eq. fl2.5p . In contrast, 
Eq. (|2.5| ) does not hold in other models with a doubly charged scalar (e.g., the Left-Right 
symmetric model in which the couplings hij are essentially arbitrary.) 

In this work we are concerned with the case of the leptonic decays of dominating, 
which is realized if hy are larger than the smallest Yukawa coupling in the SM (i.e., the 
electron Yukawa coupling, h e ~ 10™ 6 ). From Eq. fl2.5|) it follows that va ^0.1 MeV in 
this scenario, and thus the decay i? ±± — > W^W^ (which depends on v&) is negligible. 
Moreover, in the HTM the lifetime of is always short enough to ensure that it decays 
in the detectors of the Tevatron and the LHC |23[| . For very small values of i>a (e.g., 
-UA < 10 eV), the magnitude of the Yukawa couplings hij approaches unity. In this case, 
the charged scalars (which we assume to have a mass of the order of the electroweak scale) 
would induce potentially observable BR's for LFV decays such as fi — > ej, fi — > eee and 
r — > III (e.g., see pf|), whose decay rates depend on the parameters of the neutrino mass 



matrix, the absolute values of hij, and M H ±± or M H ± j9|, |10|, |17|] . Such constraints can be 



satisfied for appropriately small hij , the most severe constraint being from —> eee (which 
gives h ee h efl <^ 10 -7 ). In the parameter space 10~ 6 <^ hij <^ 10 -3 , the constraints from the 
above LFV decays are satisfied, even for values of M H ±± of the order of the electroweak 
scale. The BR of — > t^l^- depends on the six parameters of the neutrino mixing 
matrix, V , (with the dominant uncertainty arising from the unknown Majorana phases, 
(pi and 4>2), the unknown mass of the lightest neutrino (mo), the mass splittings of the 
neutrinos, and the ignorance of the neutrino mass hierarchy (normal or inverted) |J. 

Detailed studies of BR(i? ±=t g±t±) have been performed in [JTT], 0, [Uj . Notably, 
BR(F ± t± v) (in which the three flavours of neutrinos are summed over) does not depend 
on the Majorana phases, and the dominant uncertainty is from mo and the neutrino mass 
hierarchy Q. Importantly, BR(H ±:iz £ ± e ± ) ~ 100% and BR(i/ ± -> t±v) ~ 100% for a 
given lepton flavour are not possible in the HTM. Moreover, BR(.ff ±=l: — > £ =t £ =l= ) ~ 100% and 
BR(F ± -> (± v) ~ 100% are not possible even when summing over I = e,/i, although values 
as high as ~ 70% are possible in specific regions of the parameter space of the neutrino 
mass matrix. A generic model containing can have BR(if ±=l= — > ^ ± £ ± ) ~ 100% 

(summing over £ = e, fx) for appropriately chosen h^, e.g., the LR symmetric model with 
h ee , h efM , h^ lfl 3> h eT ; h^ T , h TT . In direct searches for H ±zil the derived lower limits on M H ±± 
are usually given for the extreme case of BR=100%, and this will be discussed in more 
detail in the next section. 

3. Searches for doubly charged Higgs bosons at hadron colliders 

Direct searches for have been carried out at the Fermilab Tevatron [27, 28, 29, |30| , 
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assuming the production mechanism at partonic level given by qq — > 7*, Z* — > H ++ H , 
whose cross section depends on only one unknown parameter, M H ±±. Production mecha- 
nisms which depend on the triplet VEV (qq' —?■ — > W T H ±± and fusion via VF ± *W ± * — > 
[24, p8l) are not competitive with qq — > H ++ H~~ at the energies of the Tevatron. All 



searches assume the leptonic decay mode — > £^1^, for which there are six possibilities 
(ee, /U/i, rr, e/i, er, fir). 

The CDF collaboration searched for three final states, — > e^e^, e fi , /x^/j^, 



requiring at least one pair of same-sign leptons with high invariant mass [27]. The integrated 
luminosity used was 0.24 fb _1 and the mass limits Mh±± > 133, 113, 136 GeV were obtained 
for the decay channels H ±zil — > e^e 1 * 1 , e ± // ± , /i ± ^ =t , respectively, assuming BR=100% in a 



given channel. The DO collaboration [28, |2j| searched for — > /x^/i 1 * 1 , and derived 



the mass limit M H ±± > 150 GeV p9|] using 1.1 fb 1 of integrated luminosity. The main 



difference between these searches by DO is the requirement in the most recent search [29] of 
a third /j, of opposite sign to the two same-sign /x, where the latter is assumed to originate 
from the decay of one of the pair-produced iJ ±=l= . This extra requirement suppresses 
backgrounds from 7/Z — > fJ, + ^~ and multijets, which were less than one event for the 
integrated luminosity of 0.11 fb -1 used in [28], but became non-negligible for the search in 



1 29] with 1.1 fb . The requirement of a third lepton is necessary for the future Tevatron 
searches in order to reduce the SM backgrounds. At the Tevatron, the main backgrounds 
to the three-lepton signal with £ = e,fi are from WZ and ZZ production. Two decay 
channels involving r (i? ±± — > e t , /i ± r ± ), were searched for by the CDF collaboration 
in [30], and there has been no search for — > 



As discussed in Section ||, in the HTM one has BR(^ ±± -> £ ± £ ± ) < 100% m a given 
channel, and thus the above mass limits for M H ±± (which assume BR=100%) are weakened 
when applied to the HTM. The expected number of —¥ i^l^ events scales linearly in 



BR (for searches for a single pair of same-sign leptons [27, [2^]) or quadratically in BR (for 
searches which require a third lepton or more |2£|), and so the mass limits are weakened 
accordingly. 

All the above searches at the Tevatron assume only the production mechanism qq — > 
7*, Z* ->■ H ++ H . However, the partonic process qq' — > W* — > H ±:k H T Q [||, |26j has 



a 



cross section at hadron colliders comparable to that of qq — > H ++ H~~ for M H ± ~ M H ±±, 
and thus the former will also contribute to the search for ff ±=t . In Ref. [26|, it is suggested 
that the search potential at hadron colliders can be improved by considering the following 
inclusive single ff ±=t cross section (<r#±±): 

a H ±± = a(pp,pp — > H ++ H ) + a(pp,pp — > H ++ H~) + a(pp,pp — > H H + ) (3.1) 

At the Tevatron a(pp — > H ++ H~) = a(pp — > H H + ) while at the LHC a(pp — > 
H ++ H~) > a(pp — > H H + ). These two production mechanisms have different QCD 
K factors. Explicit calculations |22|] for pp,pp — > H ++ H~~ give around K = 1.3 at the 
Tevatron and K = 1.25 at the LHC, with a dependence on M H ±±. In reality, the K 
factor for pp,pp —¥ H^H^ is expected to be very similar (but not identical) to that for 



pp,pp — > H ++ H , with some dependence on the mass splitting M H ±± —M H ±. In [26, 14] 
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the K factors are taken to be equal. We note that pp,pp — > H ++ H also receives a con- 
tribution from real photon annihilation [23|, which causes an increase in the cross section 
of around 10% at the LHC, but much less at the Tevatron. In our simulation analysis 
however we do not include this correction. 



4. Simulations of signal and backgrounds at the LHC 

Several studies have been performed to study the doubly charged Higgs in the decay channel 
— > if if = e,fi,r) at the LHC. The production mechanism qq — > 7*,Z* — > 



H ++ H— followed by decay H++H— -> i + i + i~i~ is studied in El E3, pi, p|, M p]. 



Only two among these studies also take into account the production mechanism pp — > 
w ±* H ±± H T Q jy^ followed by the decays -> if if = e,^,r) and H± -> 

if v. The LHC sensitivity to — > if if considerably extends that at the Tevatron, due 
to the increased cross sections and larger luminosities. The analysis of Ref. [|32| shows that 
can be discovered for m#±± < 800 GeV and C = 50 fb _1 , assuming BR(// ±=I: — > 
//^/i 1 * 1 ) = 100%. Importantly, all the above simulations suggest that as little as C = 1 fb _1 
is needed for the discovery of m^±± < 400 GeV if one of BR(i? ±=l: — > e^e^, e^p^, n^^) 
is large, and hence such a light would be found very quickly at the LHC. The signal 
from qq -> 7*,Z* -)■ H ++ H and decay H ++ H -> i + i + i~i~ is usually taken to be 
four leptons, which are isolated and have sufficiently large transverse energy. In Ref. |32| 
the signal is taken to be 4/x. In Ref. |3l|, where little Higgs models are considered, the 
signal is defined as Ai where £ = fj,,T (and e is not included), and five different four-lepton 
signatures are studied (one of which being 4^). In Ref. ]23| e and /x are treated as the same 
particle, and a parton-level study of the four-lepton signature is performed. In Ref. [|33| 
two signatures are defined: i) four leptons and ii) at least three leptons. It is shown that 
superior sensitivity to M H ±± is obtained for the signature of at least three leptons. 

As discussed earlier, the production mechanism pp — > W^* — > H^H^ will contribute 



to the signal for i^ ±± if three (or more) leptons are required. The simulation in Ref. [15] is 
the first study of the mechanism pp — > i/ ±± i? =F together with pp — > H ++ H , with the aim 
of improving the sensitivity to M H ±± at the LHC. In Ref. j0|] e and \x are not distinguished, 
and such an inclusive channel has the advantage of maximizing the sensitivity to M H ±± for 
a given integrated luminosity, and for the general case of BR < 100% for a given flavour 
of lepton. Both a four-lepton signature and a three-lepton signature are studied, and the 
sensitivity to M H ±± for the two signatures is compared, assuming M H ±± = M H ±. The 
three-lepton signature is defined as being exactly three leptons (3^), i.e., a fourth lepton is 
vetoed. Note that this three-lepton signature differs from that defined in the latest search 
for i? ±± at the Tevatron in which a fourth lepton is not vetoed (> 3£). 

In Ref. |LJ| it is concluded that the three-lepton signature offers considerably greater 
discovery potential for iJ ±=l= in the HTM than the signature of four leptons (note that 



the same conclusion is obtained in [33|, even without including pp — > H^H^). The 
main reason for the superior sensitivity of the three-lepton signature in [jl5| is the extra 
contribution from pp — > H^H^ (which does not contribute to the four-lepton signature). 
Although the SM background for the three-lepton signature is larger than that for the four- 
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lepton signature, in the region of high invariant mass of (relevant for Mjj±± > 200 

GeV) the backgrounds are still sufficiently small, which gives rise to superior sensitivity to 
M H ±± for the three-lepton signature. Moreover, Ref. |ll| used different sets of cuts for the 
three-lepton and four-lepton signatures. 

In Ref. [|14] , a parton-level study at the LHC was performed for the detection prospects 
of the production channel pp — > H^H^ alone, followed by the decays — > t^l^ and 
if ^ — > l^v, where both e and [i contributions are summed together in an inclusive approach 



like that in Ref. The strategy in [14] is to isolate the contribution from pp — > H^H^ 
and remove that from pp — > H ++ H , with the aim of probing the vertex H^H ^W^, 
which is present in the HTM but not in models with SU (2) singlet scalars. A cut is imposed 
on missing energy (which originates from H ± — > £^u) in order to remove the contribution 
from pp — > H ++ H . Therefore the approach of Ref. |l4| contrasts with that of (and 
our approach), where in the latter the cuts are designed to keep signal events from both 
pp — > H ++ H and pp — > H^H^ in order to optimize sensitivity to M H ±± for a given 
integrated luminosity. 
The main features of our analysis are : 

• In order to analyze the signature of > 3t as mentioned above we have to consider 



the H ±± H T W ± vertex. This vertex was not available in pythia |40| for the HTM 



model. For our analysis we have used CalcHEP [39] and incorporated this vertex in 
the model file. 

• We have included K-factors for both signal and background events. 

• We have performed a detailed realistic detector simulation using the fast detector 
simulator atlfast for both signal and background processes. 

• We have considered the > 3£ signature and compared its discovery potential with 
that for the 4£ signature at the LHC. 

Hereafter, we will refer to electrons and/or muons collectively as "leptons" (£ = e,fi). 
We will further assume the idealized case of BR(if ±± -> l±g±) = 100% and BR(H ± 
l^vg) = 100%, i.e., the decays of the charged Higgs bosons are saturated by the electronic 
and muonic modes. We do this in order to provide a simple comparison of the discovery 
potential of the two signatures under investigation. Moreover, such extreme branching 
ratios are generally used when deriving limits on M#±± from direct searches. In contrast, 
we note that representative branching ratios in the HTM were used in |l5|] , for which decay 
modes of if ±=t involving r were sizeable. Careful attention was given to secondary electrons 
and muons which originate from decays like — > followed by r — > ivv, and their 

effect on the dilepton invariant mass distribution was studied. In our analysis the decay 
modes of if ±=l= involving r are absent, and so there are no such secondary leptons. For the 
cases of BR(if ±± -> < 100% and BR(if ± -> t±v t ) < 100% (as discussed in Section | 

for the HTM) our results will need to be scaled by multiplicative factors of branching ratios. 
Moreover, for non-zero BRs of and H ± into final states which contain r leptons, the 
influence of the secondary leptons (which originate from the decay of the r leptons) on 
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the signal will need to be included. For definiteness, we take M H ±± = M H ± as the mass 
difference is fairly small for most of the parameter space in the HTM. 




Figure 1: Cross section of inclusive doubly charged Higgs bosons production (Eq. 3.1) as a function 



of M H ±±. The K-factor of the processes is taken to be 1.25 for LHC and 1.3 for Tevatron. 

The cross section for the inclusive production of doubly charged Higgs bosons at 
hadronic colliders is plotted in Figure [l] as a function of the doubly charged Higgs mass. 



In this plot we show the production cross sections as given in Eq. (3.1) for the LHC at the 
center-of-mass (CM) energy of y/s = 7, 10, 14 TeV and for the Tevatron y/s = 1.96 TeV. 
The LHC is expected to take around 1 fb _1 of integrated luminosity at y/s = 7 TeV in 
its first two years of operation. Subsequently, the machine is planned to run at the design 
energy of y/s = 14 TeV. For completeness, we also show the cross section at an interme- 
diate energy of y/s = 10 TeV because operation at this energy has been discussed. In our 
later simulations, we have assumed the CM energy of y/s = 14 TeV at the LHC. We have 
used the leading-order (LO) CTEQ6L parton distribution functions (PDF) with two-loop 
a s running, and identified both the factorization scale fif and the renormalization scale \i r 
with the partonic CM energy s. 

4.1 Framework for event generation 

The SM background processes we have considered for both the > 31 and 41 channels are: 

• ZZ with each of the Z's decaying leptonically. 

• W^Z with each of the weak gauge bosons decaying leptonically. 

• ti with t — > Wb, and VF's and b decaying (semi) leptonically. 
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• Zbb with Z and b decaying (semi)leptonically. 

• Wbb with W and b decaying (semi)leptonically. 

• Ztt with Z and t decaying (semi)leptonically. 

• Wtt with W and t decaying (semi)leptonically. 

• W^W^W^ with each of the W's decaying leptonically. 



m H ±± (GeV) 


200 


300 


400 


500 


600 


700 


800 


Events generated 


300 K 


150 K 


100 K 


100 K 


80 K 


80 K 


60 K 



Table 1: Number of signal events generated where K stands for 10 3 . These events are then scaled 
to the luminosity of £ = 10 fb . 

The setup for signal and background event generations is the following: 



Signal event generation: We have used CalcHEP v. 2. 5. 4 [39| for calculating cross 
sections. For this purpose we have implemented the relevant interaction vertices 
in the CalcHEP model files. The partonic level signal events have been generated 
using CalcHEP and then passed to PYTHIA v.6.4.21 |h| via Les Houches Event 



(LHE) interface [41] in order to include initial state radiation/final state radiation 
(ISR/FSR) effects. The number of signal events generated are given in Table |. 
These events were then scaled to the luminosity of C = 10 fb -1 . 

Background event generation: We have generated the ti, ZZ, W^Z events di- 
rectly using pythia , and the Zbb, W ± bb, Ztt, W ± tt events first using CalcHEP 
and then interfaced with PYTHIA for ISR/FSR. The number of background events 
generated are given in Table These events were scaled to the luminosity of C = 10 
fb" 1 . 



Process 


ZZ 


WZ 


WWW 


tt 


Z bb 


W bb 


Ztt 


Wtt 


Decay modes 


all 


all 


all 


SL 


SL 


SL 


SL 


SL 


Events 


1.5 M 


1.5 M 


300 K 


90 M 


1.2 M 


900 K 


150 K 


800 K 



Table 2: Number of background events generated where K stands for 10 3 and M stands for 10 6 . 
The second row corresponds to the decay modes considered, and SL indicates "semi-leptonic." 
These events are then scaled to the luminosity of C — 10 fb -1 . 



In order to make more realistic estimates of the signal and background events, we 
have further processed both of them through the fast ATLAS detector simulator atlfast 
1 43]. The resulting events have been analyzed within the ROOT framework. The detector 
simulator atlfast provides simple detector simulation and jet reconstruction using a sim- 
ple cone algorithm. It further identifies isolated leptons, photons, b and r jets, and also 
reconstructs missing energy. In our analysis for the LHC, we have assumed the CM energy 
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Process 


WZ 


ZZ 


tt 


Zbb 


Ztt 


1^66 




K- factors 


1.5 


1.35 


1.67 


2.4 


1.35 


2.57 


1.25 



Table 3: K- factors for the background and signal processes at LHC 



of 14 TeV and luminosity C = 10 fb _1 . The K-factors for signals and backgrounds that we 
have used are listed in Table [3|. 

In order to improve the search potential for the doubly charged Higgs boson at the 
LHC, we want to advocate the DO search strategy ]^9| of looking for > 3£. This is in 
contrast with the current LHC search strategy |3l]] of looking for exactly four leptons. 
Accordingly, we will present our simulation results for exactly four-lepton signature and 
> 3 lepton signature. 

4.2 Signature of four leptons 

This is the signature (4£, I = e, [i) when we have the pair production of doubly charged 
Higgs bosons via the process pp — > H^H^ followed by leptonic decays — > l^l^ . 



This signature has been simulated within the context of LHC in Refs. |32|, 15, 33]. 
We use the following preselection cuts on signal and background events [JDJ : 

• There are exactly four leptons with two for each charge sign (£ + £ + 1~ £~) in each 
event. 

• Each of the leptons has \p l T \ > 5 GeV and pseudorapidity in the range \t]\ < 2.5. 

• Amongst the four leptons, at least two of the leptons have \p¥p\ > 30 GeV. This cut 
reduces the backgrounds where the leptons originate from the semileptonic b decays 
as they tend to be less energetic. 

• Opposite-sign dilepton invariant mass cut: m f + f - > 20 GeV. This is done in order to 
suppress the backgrounds where the opposite-sign lepton pair comes from a photon. 

We impose additional cuts to further improve the signal significance: 

(a) The Z window cut. The invariant mass of opposite-sign dileptons is required to be 
sufficiently far from the Z mass: \rri£±£^ — Mz\ > 10 GeV. This removes events where 
the leptons come from the Z decay. 

(b) The Ht cut. One can also use the total transverse energy (H?) as a parameter to 
distinguish signals from backgrounds. The total transverse energy is defined as 

Ht = Y, \Pr\ ■ 

The Ht distribution of signals tends to peak around the heavy particle mass, and 
hence this cut can be used to suppress the SM background processes that involve 
relatively light particles. We will show the results for Ht > 300 and 500 GeV. We 



note that this Ht cut was also used in the study of the 4£ signature in Ref. [31], but 
it was not used in the study of Ref. [0 . 
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Backgrounds 


Signal (M H ±±) 




WZ 


ZZ 


tt 


Zbb 


Ztt 


Wtt 


200 GeV 


600 GeV 


Pre-selection 


0.2 


130.5 


1.3 


0.2 


122.6 


0.1 


400.1 


4.2 


\me+e- — m z\ > 10 GeV 


0.1 


2.1 


0.3 





2.1 


0.1 


330.6 


4.1 


H T > 300 GeV 





0.4 








1.2 





327.9 


4.1 


H T > 500 GeV 





0.1 








0.3 





222.9 


4.1 


S 














48.7 


3.7 



Table 4: Background and signal events surviving the cuts for exactly 4-lepton final states. For 
these numbers we have taken C = 10 ftr 1 and = 14 TeV. 



The effects of the cuts on signal and background events for two indicative masses of 
the doubly charged Higgs boson are given in Table ||. After the pre-selection cuts, the 
dominant backgrounds are ZZ and Ztt, which are then significantly reduced by the Z 
window cut. For large Mjj±±, the Z window and Ht cuts only reduce the number of signal 
events by a negligible amount and hence the significance of the signal is improved. 



4.3 Signature of at least three leptons 







D/f 




□ zz 




rjzu 




1 \Zbb 




m H* = 


200 GeV 


m H" = 


400 GeV 


m H* = 


600 GeV 




300 400 

U T (GeV) 



Figure 2: fix distribution for doubly charged Higgs 
production (signal) and SM backgrounds after imposi- 
tion of the pre-selection cuts. We have taken the CM 
energy of y/s = 14 TeV and luminosity C = 10 fb _1 . 



The study of the signature with ex- 
actly three leptons (3^) was done in 



Ref. [15], with the motivation of com- 
paring the detection prospects of three 
distinct types of seesaw-based models 
of neutrino mass generation (one of which 
being the HTM). It was acknowledged 
that such a signature is not necessarily 
the one which optimizes the discovery 
potential in a given model. In this pa- 
per, we are concerned with optimizing 
the sensitivity to M H ±± and thus we 
consider the signature of > 2>l as done 



in the DO search [29]. As emphasized 
earlier, both the pair production and 
single production of the doubly charged 
Higgs boson contribute to the signature 
and hence increase the signal events. 
For this analysis, we have used cuts 
similar to those in Section 4.2. The 



only difference is that instead of singling out events with exactly four leptons, here we 
select events that have at least three leptons, of which two have the same sign. We would 



like to note that Ref. [15| used a slightly different set of pre-selection cuts for studying 
the exactly three-lepton signature. For the exactly three-lepton signature, they chose events 
that had at least two same-sign leptons with px > 30 GeV. Also, the pre-selection cuts used 
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Backgrounds 


Signal (M H ±±) 


Pn-t-e II 

ours -\y 


WZ 


WWW 


ZZ 


tt 


Zbb 


Ztt 


Wtt 


200 


600 


Pre-selection 


591.7 


3.5 


203.6 


159.9 


57.7 


212.5 


9.7 


1570.4 


17.6 


\me+e- ~ m z\ > 10 GeV 


50.9 


2.7 


12.1 


113.2 


0.9 


33.4 


7.4 


1397.8 


17.3 


H T > 300 GeV 


7.5 


1.1 


1.6 


8.9 





17 


3.4 


1351.1 


17.3 


H T > 500 GeV 


1.7 


0.3 


0.4 


0.9 





3.2 


0.6 


796.2 


17.3 


S 
















77.4 


5 



Table 5: Background and Signal events surviving the cuts for at least 3 leptons in the final state. 
We have taken £ = 10 ftr 1 and ^ = 14 TcV. 



in Ref. [15] were slightly different for exactly four-lepton and exactly threedepton signa- 
tures, whereas in order to make a comparison between the fourdepton and > threedepton 
signatures we are using the same set of cuts for both. 

The effects of the cuts on signal and background events for the > 3£ signature are 
given in Table ||. It is evident that the discovery reach of the current signature is better 
than the 4£ signature. 

Also shown in both Tables |] and |B| are the significance. As signal and background 
events become fewer after the cuts, it is necessary to employ Poisson statistics to estimate 
the significance of the signal. We use the significance estimator [44| 



2<jn ln(l + ^) -s}, (4.2) 



b 

where b is the expected number of background events and no is the number of observed 
events. Accordingly, the number of signal events is s = tiq — b. This estimator is based 
on a logdikelihood ratio, and follows very closely the Poisson significance. In the limit of 
s/b 1, it reduces to the simple estimator s/y/b. 

The distributions presented in Figures ^, |3|, and || are for the luminosity of C = 10 fb _1 
and CM energy of 14 TeV at the LHC. In Figure we show the $t distribution for both 
signal and background events. The background events concentrate at lower whereas 
the signal events extend to the higher region as they contain more energetic neutrinos 
coming from the singly charged Higgs bosons. In Figure ||, we show two dilepton invariant 
mass distributions. As can be seen in the opposite-sign dilepton invariant mass distribution 
(left panel), the event distribution tends to peak around the Z mass as the lepton pairs 
mostly originate from the decay of the Z boson. These backgrounds can be readily reduced 
by imposing the Z window cut on the opposite-sign dilepton invariant mass. On the other 
hand, the same-sign dilepton invariant mass distribution peaks around the doubly charged 
Higgs mass whereas the SM processes form a continuous background, as shown in the right 
panel of Figure |3| 

Finally we show the total transverse energy distribution in Figure f§ The H T distri- 
bution peaks around the total mass of the heavy particles produced in the hard process. 
Hence this variable serves as a very useful discriminator between the signal and SM back- 
grounds, especially for relatively heavy doubly charged Higgs masses. We emphasize that 
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Figure 3: (a) Left panel: opposite-sign dilepton invariant mass distribution, (b) Right panel: 
same-sign dilepton invariant mass distribution, for both signal and background events. 




Figure 4: Total transverse energy (Ht) distribution for signal and SM background events. 



we have performed the first study of the dependence of the three-lepton signature on the 
parameter Ht- It is evident that a clear signal for H^ can be obtained by using a cut on 
Ht to reduce backgrounds, even before plotting the dilepton invariant mass distribution 
(the latter providing information on M H ±±). 
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Discovery potential 




200 300 400 500 600 700 

m H±± (GeV) 



Figure 5: Luminosity required for a 5cr discovery of the doubly charged Higgs boson as a function 
of its mass at the LHC. The two curves correspond to the exactly 4-lepton signature and the > 3£ 
signature, respectively. 

5. Discussions and Summary 

In this paper we have analyzed the leptonic signatures of the production of a doubly charged 
Higgs boson at LHC. For this purpose we have used CalcHEP to generate the signal events, 
and then interfaced it with PYTHIA . For more realistic estimates of signal and background 
events, we have used fast ATLAS detector simulator atlfast . We have also included 
relevant K-factors for both signal and backgrounds in our analysis. 

Using the significance estimator, we have estimated the LHC discovery potential for the 
doubly charged Higgs boson in the 4-lepton mode from the pair production only and the > 3 
lepton mode from the inclusive production, i.e., both pp — > H ++ H and pp — > H^H^. 
We have performed the first simulation of the > 3 lepton channel. Our result is shown 
in Figure || where the luminosity required to make a 5a discovery of the doubly charged 
Higgs boson is plotted as a function of its mass. It is clear that the discovery potential for 
i? ±± at the LHC through the > 3£ mode is significantly better than the 4-lepton mode. 
For an integrated luminosity of 10 fb , for example, one detector at LHC alone can reach 
~ 600 GeV for the former and ~ 550 GeV for the latter. 

In our analysis of multi-lepton signatures we have not considered the QCD background 
where a jet is misidentified as a lepton. Although the probability of misidentification is quite 
small (3^], the QCD production cross-sections are many orders of magnitude larger than 
the multi-lepton (3^ and Ai) cross-sections and hence can contribute to the backgrounds. 
Importantly, the QCD background for the four-lepton signal will be smaller than that for 
the three-lepton signal. Hence the inclusion of the QCD background would introduce a 
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systematic error into the estimates given in Figure [|, which could alter the significance of 
our results. In an actual experiment, the QCD background can be estimated from the data. 
For a more realistic study, one should do a full detector simulation, optimize the cuts, and 
take into account the statistical and systematic uncertainties. We hope that this study 
will motivate our experimental colleagues at hadronic colliders to update their analyses by 
considering the following points: 

• Tevatron: Include the process pp — > H^H^ in the analysis when searching for the 
> 3£ signature. 

• LHC: Search for the > 31 signature to increase the LHC discovery reach of the 
doubly charged Higgs boson. 

Acknowledgments 

C. C. would like to thank Neil Christensen for valuable help with the CalcHEP pack- 
age. N. G. would like to thank Michael Klasen and Arun Kumar Nayak for their com- 
ments/suggestions, and acknowledge the support from NCTS, Taiwan for his visit to NCU, 
where this work was initiated. This work of C. C. was supported in part by the Grant 
No. NSC 97-2112-M-008-002-MY3 and NCTS. A. G. A was supported by the "National 
Central University Plan to Develop First-class Universities and Top-level Research Cen- 
ters", and by a Marie Curie Incoming International Fellowship, FP7-PEOPLE-2009-IIF, 
Contract No. 252263. N. G. was supported by a grant from University Grants Commission 
(UGC), India under project number 38-58/2009 (SR). 

References 

[1] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1562 (1998). 

[2] Y. Kuno and Y. Okada, Rev. Mod. Phys. 73, 151 (2001); M. Raidal et ai, Eur. Phys. J. C 
57, 13 (2008). 

[3] W. Konetschny and W. Kummer, Phys. Lett. B 70, 433 (1977). 

[4] R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

[5] M. Magg and C. Wetterich, Phys. Lett. B 94, 61 (1980); G. Lazarides, Q. Shah and 
C. Wetterich, Nucl. Phys. B 181, 287 (1981). 

[6] J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980). 

[7] T. P. Cheng and L. F. Li, Phys. Rev. D 22, 2860 (1980). 

[8] E. Ma and U. Sarkar, Phys. Rev. Lett. 80, 5716 (1998); E. Ma, M. Raidal and U. Sarkar, 
Phys. Rev. Lett. 85, 3769 (2000); E. Ma, M. Raidal and U. Sarkar, Nucl. Phys. B 615, 313 
(2001). 

[9] E. J. Chun, K. Y. Lee and S. C. Park, Phys. Lett. B 566, 142 (2003). 
[10] M. Kakizaki, Y. Ogura and F. Shima, Phys. Lett. B 566, 210 (2003). 
[11] J. Garayoa and T. Schwetz, JHEP 0803, 009 (2008). 



- 15 - 



A. G. Akeroyd, M. Aoki and H. Sugiyama, Phys. Rev. D 77, 075010 (2008). 

M. Kadastik, M. Raidal and L. Rebane, Phys. Rev. D 77, 115023 (2008). 

P. Fileviez Perez, T. Han, G. y. Huang, T. Li and K. Wang, Phys. Rev. D 78, 015018 (2008). 

F. del Aguila and J. A. Aguilar-Saavedra, Nucl. Phys. B 813, 22 (2009). 

S. T. Pctcov, H. Sugiyama and Y. Takanishi, Phys. Rev. D 80, 015005 (2009). 
A. G. Akeroyd, M. Aoki and H. Sugiyama, Phys. Rev. D 79, 113010 (2009). 

A. G. Akeroyd and C. W. Chiang, Phys. Rev. D 80, 113010 (2009). 

T. Fukuyama, H. Sugiyama and K. Tsumura, JHEP 1003, 044 (2010). 

V. D. Bargcr, H. Baer, W. Y. Keung and R. J. N. Phillips, Phys. Rev. D 26, 218 (1982). 

J. F. Gunion, J. Grifols, A. Mendez, B. Kayser and F. I. Olness, Phys. Rev. D 40, 1546 
(1989); J. F. Gunion, C. Loomis and K. T. Pitts, eConf C960625, LTH096 (1996) 
[arXiv:hep-ph/9610237]. 

M. Muhlleitncr and M. Spira, Phys. Rev. D 68, 117701 (2003). 

T. Han, B. Mukhopadhyaya, Z. Si and K. Wang, Phys. Rev. D 76, 075013 (2007). 

K. Huitu, J. Maalampi, A. Pietila and M. Raidal, Nucl. Phys. B 487, 27 (1997); J. Maalampi 
and N. Romancnko, Phys. Lett. B 532, 202 (2002). 

B. Dion, T. Gregoire, D. London, L. Marleau and H. Nadeau, Phys. Rev. D 59, 075006 
(1999). 

A. G. Akeroyd and M. Aoki, Phys. Rev. D 72, 035011 (2005). 

D. E. Acosta et al. [CDF Collaboration], Phys. Rev. Lett. 93, 221802 (2004). 

V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 93, 141801 (2004). 

V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 101, 071803 (2008). 

T. Aaltoncn et al. [The CDF Collaboration], Phys. Rev. Lett. 101, 121801 (2008). 

A. Hcktor, M. Kadastik, M. Muntcl, M. Raidal and L. Rebane, Nucl. Phys. B 787, 198 
(2007). 

T. Rommcrskirchcn and T. Hebbeker, J. Phys. G 34, N47 (2007). 

G. Azuelos, K. Benslama and J. Ferland, J. Phys. G 32, 73 (2006). 

B. Pontecorvo, Sov. Phys. JETP 7, 172 (1958) [Zh. Eksp. Teor. Fiz. 34, 247 (1957)]; Z. Maki, 
M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

S. M. Bilenky, J. Hosek and S. T. Petcov, Phys. Lett. B 94, 495 (1980); M. Doi, T. Kotani, 

H. Nishiura, K. Okuda and E. Takasugi, Phys. Lett. B 102, 323 (1981). 

F. Cuypers and S. Davidson, Eur. Phys. J. C 2, 503 (1998). 

T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 10, 113011 (2008). 

R. Vega and D. A. Dicus, Nucl. Phys. B 329, 533 (1990). 

A. Pukhov, arXiv:hep-ph/0412191. 

T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) [arXiv:hep-ph/0603175]. 



- 16 - 



[41] J. Alwall et at, arXiv:0712.3311 [hep-ph]; E. Boos et ai, arXiv:hep-ph/0109068 ; J. Alwall et 
at, Comput. Phys. Commun. 176, 300 (2007) [arXiv:hep-ph/0609017]. 

[42] F. Maltoni and T. Stclzer, JHEP 0302, 027 (2003) [arXiv:hep-ph/0208156]. 

[43] E. Richter-Was et.ai. , ATLFAST 2.2: A fast simulation package for ATLAS, 
ATL-PHYS-98-131. 

[44] G. L. Bayatian et al. [CMS Collaboration], J. Phys. G 34, 995 (2007). 



-17- 



